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ABSTRACT 


This study examines enlisted attrition behavior for the 
U.S. Coast Guard and develops a model that projects attrition 
figures. Survival analysis techniques are used to analyze the 
empirical attrition behaviors associated with an individuals 
sex, race, marital status, and military occupational skill 
(MOS). In this study males tend to have higher survival 
probabilities than females, non-caucasians higher’ than 
caucasians, and married persons higher than those not married. 
Aviation MOSs have the highest survival probabilities and 
technical MOSs have the lowest. 

Modelled survivor functions are developed for two 
paygrades because each contain small personnel inventories. 
These modelled survivor functions do not fit the data as well 
as desired but are nonetheless used pending the development of 
sharper alternatives. 

Finally a counting model based on the Binomial 
Distribution is developed that projects monthly enlisted 


attrition figures. 
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I. INTRODUCTION 


A. PROBLEM STATEMENT 


The United States Coast Guard (CG) is a military service 
operating within the Department of Transportation. Its main 
missions are aids to navigation, boating safety, defense 
operations, maritime law enforcement, environmental response, 
marine inspection, ice operations, marine licensing, marine 
science, search and rescue, port security and safety, and 
waterways management. The CG maintains an active duty 
hierarchical personnel structure with approximately 6,000 
officers and 30,000 enlisted members. The CG Personnel 
Workforce Planning (PWP) office located in Washington, D.C. 
has the overall responsibility for projecting the active duty 
enlisted personnel needs of the CG. Their responsibilities 
include the calculating of personnel stocks, promotion 
requirements, attritions, and recruitment needs. The current 
method for predicting monthly enlisted attrition by paygrade 
is to calculate a running average or mean of past monthly 
attrition and use this mean as a forecast of future attrition. 
This forecasting method is crude, ignores many issues, and is 


not accurate enough to meet the needs of the CG. 


Why does a person separate from the CG? What are the 
characteristics associated with this person? Are these 
characteristics the same for each paygrade? Previous CG 
studies have focussed on the first question but have not 
looked at attrition data analytically. The first question 
will not be addressed. This thesis will look at the next two 
questions and will propose a method, based on a survival 
analysis of CG enlisted attrition data, for calculating 
monthly personnel attrition. This study is a first attempt at 
examining CG enlisted attrition using data analysis tools. It 
is designed to provide a procedure for studying attrition, and 
for developing a basic model that calculates monthly attrition 


which can be expanded by future studies. 


B. COAST GUARD ENLISTED PERSONNEL BACKGROUND INFORMATION 
The Coast Guard enlisted personnel structure is composed 
of 9 paygrades, E-1 to E-9, lowest to highest respectively. 
(There is an E-10 position filled by one person but this 
thesis does not address this paygrade.) Paygrade E-1 is 
composed of persons attending CG basic training school and 
paygrade E-2 is composed of persons who have completed basic 
training school, have been assigned to active duty commands, 
and are preparing to attend a CG specialty school (also called 
A-school). Paygrades E-3 through E-9 are divided into 24 
active duty military occupational skills (MOS) or _ sub- 


specialties listed below: 


Aviation Machinist’s Mate 
Aviation Electrician’s Mate 
Aviation Structural Mechanic 
Aviation Survivalman 
Aviation Electronic Technician 
Boatswain’s Mate 

Damage Controlman 
Electrician’s Mate 
Electronics Technician 

Fire Control Technician 
Gunner’s Mate 

Health Services Technician 
Machinery Technician 
Musician 

Marine Science Technician 
Public Affairs Specialist 
Quartermaster 

Radarman 

Radioman 

Storekeeper 

Subsistence Specialist 
Sonar Technician 

Telephone Technician 


Yeoman 


(AE) 
(AM) 
(ASM) 
(AT) 
(BM) 
(DC) 
(EM) 
(ET) 
(ED) 
(GM) 


(HS) 


(MU) 
(MST) 
(PA) 
(QM) 
(RD) 
(RM) 
(SK) 
(Ss) 
(ST) 


Cr) 


An enlisted man (EM) receives an MOS after either 
completing a specialty school or by completing an intensive 
on-the-job training program. Both methods are usually 
completed while the EM is in the E-3 paygrade. An E-3 who has 
not obtained an MOS is called nonrated; E-3’s through E-9’s 
hold only one MOS at a time and usually keep that MOS 
throughout their careers. 

Promotions to the next highest paygrade are performed by 
MOS and are dictated by the needs of the CG. Promotions to 
paygrades E-5 through E-9 are vacancy driven and promotions to 
paygrades E-3 and E-4 occur on qualification for advancement. 
Promotion to paygrade E-2 occurs upon completion of basic 
training. 

Enlisted personnel sign service obligation contracts 
(enlistments) that require the person to serve in the CG for 
a set number of years. Enlistments usually last four years. 
(There have been two year enlistments but these have been 
terminated. ) Upon completion of an enlistment and upon 
approval by the CG, a person may sign a new contract 
(reenlist), or Separate from the CG. Enlisted attrition 


generally occurs because of: 


- retirement - minimum 20 years active duty service, 


- non-reenlistments - an EM chooses not to reenlist or the 


CG chooses not to reenlist the person, 


- administrative reasons - an EM may depart the CG prior to 
the end of his/her enlistment "for the convenience of the 


government", 


- death or disability, 


- selection to an officer program. 


For the purposes of this study the term separation refers to 


a person leaving the CG for any of the above reasons. 


or. PROCEDURE 

This thesis examines, by paygrade, individual personnel 
characteristics associated with attrition and develops a model 
that calculates monthly attrition. Survival analysis 
techniques are used to study behaviors associated with these 
attrition characteristics. Survival analysis is a data 
analysis technique that is useful in examining failure time 
data. Many of the visible applications have been in medical 
research where a typical "failure" is the death of a patient 
who had been subjected to some type of treatment. Obviously 
it would not be practical to wait for all the patients to 
"fail" before analyzing the data. It is more convenient to 
end the experiment at some specified time and do an analysis 
of all the data. The observables that did not "fail" by the 


end of the data collection period are called censored 


observations. They must be treated differently since their 
failure times are yet to occur and are hence unknown. 
Survival analysis allows for censored failure times and 
provides techniques for including them in the analysis. The 
primary tool used in the survival analysis will be the 
estimation of survivor functions or curves. These curves 
provide probabilities that a person survives up to and beyond 
a set time, t. Their estimators are often called product 
limit estimators because of the technique used. Survival 
probabilities are the compliment of attrition rates and are 
needed to calculate expected monthly attrition figures. 

The personnel data are cross classified by paygrade and 
each paygrade is partitioned by certain individual personnel 
characteristics. Paygrades E-1 and E-2 will be partitioned by 
the characteristics of sex, marital status, and race. 
Paygrade E-3 will be partitioned by these same characteristics 
and by whether the person had obtained an MOS or was nonrated. 
Paygrades E-4 through E-9 will be partitioned by the same 
characteristics as E-1s and E-2s and by CG MOS. MOSs will be 
grouped according to categories listed in Table 1 on 


page 7. 


TABLE 1. MOS GROUPINGS 


oe a ————— ————————————————————— 


AD, AE, AM, ASM, AT 
Administrative 


r 
| 


Technical Eee yer TyeMsStTs STy err 


Skilled HS, SS, MU, PA 
CG Professional BM, DC, GM, MK, QM, RD, RM 





MOSs are grouped for two reasons. MOSs within the same 
group have similarities that should cause them to exhibit 
common attrition behavior and hence have common survival 
probability values. The MOSs within each category all have 
several common attributes; specifically, the type of training 
needed to acquire the MOS, the type of duty or jobs associated 
with those MOSs, and the possible civilian employment 
opportunities available for persons with these MOSs. MOSs are 
also grouped to increase the personnel inventory within each 
paygrade required for a survival analysis. Chapter II 
contains an explanation on these MOS groups and examines the 
data used in this study. 

For most MOSs stocks are large and empirical survival 
probabilities can be used. Modelled survival probabilities 
are needed if personnel stocks are smaller than preferred. 
Small personnel stocks are subject to large data variability; 


models are developed to provide stability. Paygrades E-8 and 


E-9 could be considered too small to use empirical survival 
probabilities, thus regression models are examined in Chapter 
Iv for these paygrades. 

A model that projects an expected number (and variance) of 
monthly attritions will be developed and demonstrated in 
Chapter V. This model is based on principles of the Binomial 
distribution. Paygrades will be partitioned several times to 
produce multiple binomial experiments from which expected 
values and variances can be calculated. The product limit 
estimates generated by the survival analysis will correspond 
to the probabilities of success associated with a binomial 
experiment. 

Finally conclusions and recommendations are contained in 
Chapter VI. Since this is the first time in recent years that 
enlisted attrition has been viewed using operations analysis 
tools, most of the conclusions should be unique. It is 
expected that additional studies utilizing more refined data 
will be needed to better understand, analyze, and predict CG 


enlisted attrition. 


II. DATA OVERVIEW 


A. POPULATION 

Data for this thesis were provided by the CG Pay and 
Personnel Center located in Topeka, Kansas. The data were 
provided in two sets. The first set contained individual 
personnel records for all persons who separated during FY 91 
(1 November 1990 - 31 October 1991). The second set contained 
individual personnel records for all enlisted persons who were 
on active duty the last day of FY 91 (31 October 1991). This 
second data set is in a sense a "snapshot" of the GC EM 
GCOMmpeosztlon onj31 October 1991. In both data sets each 
personnel record contained the persons sex, MOS (if 
applicable), rank, marital status, ethnic origin, date of 
current enlistment, and date of separation (if applicable). 
Each persons length of service (in months) can be calculated 
from the data entries date of current enlistment and date of 
separation. It should be noted that the entry, date of 
current enlistment, did not necessarily correspond to when the 
person first joined the CG. This affects the calculation of 
a persons time in service. In certain cases the CG allows an 
EM to separate from the CG and then allows the person to 
rejoin the CG in the future at his or her past paygrade. Such 


persons’ past service time is not ignored. The person 


maintains all accumulated service time for pay and retirement 
purposes. This condition was not reflected in the data 
provided for this study. There are numerous records of 
senior EMs having relatively little time in service when in 
all likelihood they have quite a bit of service time. This 
Situation 1s most prominent in the survival analysis results 
of Chapter III. Appendix A contains procedures and Fortran 


programs needed to access and use the data. 


B. LIMITATIONS 

Analyzing only one years amount of personnel data presents 
two problems. The first and foremost is that the data lack 
temporal stability. More than one years worth of data is 
required to provide stability. Attrition study conclusions 
and modelling results based on unstable data are not as 
reliable as if they were based on time tested data. Secondly, 
the data may not have been as reliable as preferred because 
data inputs could not be verified. The FY 91 data were placed 
in storage by the Coast Guard PPC and were not accessible 
during the study. Recommendations contained in Chapter VI 


address these shortcomings further. 


C. MOS GROUPS 
As stated in Chapter I MOS groups were formed because of 
common attributes associated within each group. These 


attributes are listed in this section. 
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1. Aviation 
Individuals with aviation MOSs have common A-school 
characteristics and common duty stations. (The CG has three 
basic types of duty stations: ashore, afloat, and air 
Stations.) The aviation A-schools (specialty school) waiting 
lists are relatively long compared to the other MOS groups. 
This lengthens a persons service time in paygrade E-3, thus 
increasing the amount of time needed for promotion to E-4. 
Unlike many MOSs, individuals within this group are limited in 
their types of duty stations. They spend most of their entire 
career at CG air stations. 
2. Administrative 
These MOSs are again grouped because of their A-school 
characteristics, choice of duty locations, and possible 
Civilian employment opportunities. The administrative A- 
schools have short waiting lists and short school durations. 
Promotion to E-4 in this MOS group is relatively fast compared 
to other groups. Administrative persons also can serve in 
almost any CG duty station or command. Since these persons 
skills are mainly clerical in nature it is felt that their 
Civilian employment opportunities are limited. 
3. Technical 
Persons with technical MOSs share the same type of 
characteristics as the administrative MOS group, but with two 


opposite attributes. Their A-schools have long waiting lists, 


ee 


are long in duration, and have low graduation rates (as 
compared to the other MOS group A-schools). Like the 
administrative group, persons in this group can serve in 
almost any CG command. However, unlike the above group, 
technically grouped persons probably have more civilian job 
opportunities. 
4. Skilled 

This group is characterized by their scope of possible 
Civilian employment opportunities. Persons within this group 
probably have marketable skills, possibly not as marketable as 
the technical group but more marketable than any other rate 
Girenor 

5. CG Professional 

This group is characterized by MOS obtainment and 
employment opportunities. As discussed in Chapter I, MOSs are 
obtained by graduating from an A-school or by completing an 
on-the job training program. This group is the only MOS group 
that incorporates on-the-job training programs for a number of 
their MOSs. MOSs are also designed for CG specific functions 
(ie. shipboard machinists, CG weapons technicians, Shipboard 
Navigators, etc.). These skills have limited value for 


Civilian employment. 
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D. SUMMARY OF PERSONNEL STOCKS 
Table 2 provides a summary of personnel = stocks 
partitioned by the rank, sex, race, and marital status for all 


CG EMs who were on active duty as of 31 October 1991. 


_ TABLE Dass ACTIVE DUTY PERSONNEL STO‘ DUTY PERSONNEL STOCKS FY 91 


fe tee ie 

137 : 
es ee eee 
e-3 [3651 |3358__|293 [3031 _|620 [1256 | 2395 
| s-4 |706a | 6364 | 700 | 6043 | 1021 | 3641 | 3423 
s-¢ |ss0a_|s4aio_|4oa | 4930 | 1074 | 4278 [1626 
s-6 |s748 [5427 |321  |4e78 [870 |4s05 |o43 
ps7 {arse ness leo jase [20 _{2ase _1az8 __ 
51752 
eum PR 














Note that race is divided into two categories: caucasian 
and non-caucasian. There are relatively few (as compared to 
the total EM population) minorities inthe CG. All minorities 
were grouped to increase the personnel stocks within the race 
category. As can be seen the CG is a male, caucasian 
dominated service. Table 3 on page 14 provides the same 


statistics for Coast Guardsmen who separated during FY 91. 


ae 


TABLE i 3. PERSONNEL STOCKS FOR "OR PERSONS § SEPARATING DURING FY 91 


lbmade Cauc. Married 

jes ba _far haere w foe 
sz tee [ee hen ol ee ee = a 
je31_[731_ [aso [72 |ar9 [276 sss 
jizae_|1094 [254 2079 |169 [602 fase 
ja3a__|308 [40 [402 [36 _—[2er_—as7_ 
rs-s | 205) af-se9 ise apoio is omegaaaa 
fe-7_|315 [305 [20 |27 [aa [277 3a 
fees pao 97 he sha ee a 
fe-9 [ae [a6 |o Jas [2 ~~ |39_ 



















Approximately the same proportions of males to females, 
caucasians to non-caucasians, and married persons to persons 
not married separated from the CG as did those who did not 
separate during FY 91. 

Table 4 on page 15 provides on a summary of personnel 


stocks for FY 91 active duty personnel by MOS group. 
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TABLE 4. FY 91 ACTIVE DUTY PERSONNEL STOCKS BY MOS GROUP 


paygrade Skilled 


ts Tal 





The numbers in paygrade E-3 are small because not all E-3s 
have an MOS. Overall approximately 48% of all personnel 
belong to the CG professional MOS group. About 14% belong in 


the aviation group; 15% belong in the administrative group; 


14% belong in the technical group; and 9% belong in the 
skilled rate group. Percentages within each paygrade are 
roughly the same. Table 5 on page 16 provides the same 


Statistics for personnel who separated during FY 91. 


LS 


TABLE 5. FY 91 SEPARATION PERSONNEL _STOCKS BY MOS GROUP 


Pestana Skilled 





hs ee ee 
a 


E. DATA PARTITIONING 


The data was partitioned by rank to form 9 paygrades. 
Each paygrade was further divided by the personnel 
characteristics sex, race, marital status, and MOS groups 
listed in Table 1 on page 7. 

1. Paygrades E-1 and E-2 

These paygrades were divided by the sex, race, and the 
marital status of the individual. Marital status indicated 
whether or not a person was currently married. 

2. Paygrade E-3 

In addition to the above, this paygrade was 
partitioned by the obtainment of an MOS. The person either 
had obtained one of the MOSs listed in Chapter I or had not 


(nonrated). 
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3. Paygrades E-4 through E-6 
These paygrades were also partitioned by sex, race, 
and marital status. Additionally they were partitioned by MOS 
groups. 
4. Paygrades E-7 through E-9 
These paygrades were divided by race, marital status, 
and MOS group. Sex was not examined because of the scarcity 


of females in these paygrades. 


ay 


III. SURVIVAL ANALYSIS 


A. INTRODUCTION 

As stated in Chapter I, survival analysis provides a 
method for analyzing censored failure time data. Failure 
times in this study refer to the number of months a person 
serves on active duty prior to separation. Failure times are 
assumed to be discrete and are rounded up to the next whole 
month figure. For example, if a person separates after 31 
months, six days of service his or her failure time is 32 
months. This chapter reviews the theory of survival analysis 
and graphically compares survivor functions of individual 
paygrades that have been partitioned by single personnel 
characteristics. Appendix C contains the SAS LIFETEST 
procedure and Fortran code required to estimate the survivor 


FuNnctE VTonce 


B. COMPUTATION OF SURVIVOR FUNCTIONS 

Let T be a nonnegative random variable representing the 
separation time of an EM from a homogeneous population. The 
survivor function is defined to be the probability that T is 


at least as great as some time, t, 
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Seer (Tees 0 <°b-<3. Cy) 
From a simple random sample, this can be estimated by 


ee oe number of observations with failure time < Cy 


Be 
total number of observations 


If T is a discrete random variable taking values x, < x, ... < 


x, with associated probability function f(x) = Pr(T=x), then 
7 
i=l 


where H is the Heaviside function, H(x)={0 if x<0; 1 if x20} 
Reread pp. 5-7). 

The presence of censored data, however, weakens the 
estimation of the survivor function in the right tail. All 
that can be safely assumed about an observation censored at 
time t is that the unobserved failure time is greater than t. 
This condition is relieved by using an alternate formulation 
of the survivor function using hazard rates. For discrete 
random variables, the hazard rate, X,, is the probability that 


the observation will immediately fail at time t given is has 


no 


Survived to time t. In the discrete case this can be 


estimated by 


s - number of observations that fail att 
e 


a noe mE EEE PT RE ee 
number of observations with failure time> t ( 


Using hazard rates [Ref. 3], the survivor function may be 


expressed as 


s(t)= [J] (1-A,). (Siyb)) 


tix 


where A; is the hazard rate at x for 1-77 ° 223). Bstimatinge 
would allow an estimation of the survivor function, S(t). 
Chapter 2 of Reference 1 contains a detailed explanation of 


the computation of the estimated survivor function. 


C. SURVIVAL ANALYSIS RESULTS 

In all paygrades the marital status characteristic and sex 
characteristic behaviors were identical; married persons 
appeared to have a higher probability of survival than those 
who were not married, and males tended to have higher survival 
probabilities than did females. (See Figures 1 and 2 on page 
40.) Non-caucasians in paygrades E-3 through E-9 had higher 
Survival probabilities than did caucasians. Results for other 


characteristics are given below. 


20 


ee Sn 1 
Caucasians and non-caucasians have about the same 
Survivor functions until about two months. After two months 
Caucasians tend to have higher survival probabilities. (See 
Figure 3 on page 41.) 
2. E-2 


Caucasians had higher survival probabilities than did 


non-caucasians up to about 20 months. (See Figure 4 on page 
41.) After that non-caucaSians had higher’ survivor 
probabilities. Also note the significant drop in 


probabilities at 24 months. The CG for a brief time allowed 
two year enlistments. This drop corresponds to enlisted 
persons separating at the end of a 2 year enlistment. 
a. E-3 
Individuals possessing an MOS had higher survival 
probabilities than did nonrates after 22 months. (See Figure 
5 on page 42.) 
4. E-4 and E-5 
The aviation MOS group in both paygrades had the 
highest survival probabilities and the technical group had the 
lowest. The remaining three groups survival functions fell in 
between the aviation and technical survivor functions. (See 


Figure 6 on page 42.) 


Zo 


5. E-6 
All MOS groups in this paygrade were very similar as 
shown in Figure 6 on page 42. All had very high probabilities 
of survival and all were very similar in shape. Again the 
aviation group had the highest survival probabilities. 
6. E-7, E-8, and E-9 
The survivor functions in all three paygrades appeared 
fairly flat with slightly decreasing slopes until about 220 
months. (See Figure 7 on page 43.) After 220 months all 
Survivor functions decreased at about the same rate. All 
survivor functions within these paygrades are characterized by 
high survival probabilities with the aviation group having the 


highest probabilities of survival. 


ae 


IV. MODELLED SURVIVOR FUNCTIONS 


A. INTRODUCTION 

Regression models are typically used to model survival 
experiences from homogeneous populations [Ref. 1]. These 
models allow for failure times to be dependent on observed 
explanatory variables. These explanatory variables can be 
either qualitative or quantitative. Two common regression 
models are the exponential regression model and the weibull 
regression model. Weibull regression models will be 
constructed and analyzed for paygrades E-8 and E-9. The 
primary reference for this chapter is Chapter 2 and Chapter 4 


of Reference 1. 


B. WEIBULL REGRESSION MODEL 
1. Weibull Distribution 
The Weibull distribution is a generalization of the 
exponential distribution that allows for a power dependence, 
p, of the hazard of time. The Weibull distribution is in fact 
an exponential distribution when p=1. The probability density 


function (pdf) of the Weibull is 


MCE SND (AES a7") "for £20 . (ea) 


The hazard function is 


ZS 


KENDAL) Ee (An) 


and the survivor function is 


StQicoGto om. (4.3) 


A convenient check to verify that data fits the 
Weibull distribution is to plot the log[-log S(t)] versus log 
t(S(t) is the empirical estimate of the survivor function) . 
If the data can be fitted by the Weibull distribution, the 
plot should provide a relatively straight line with the slope 
of the line roughly estimating the value of p. 

2. Log-Linear Regression Model 

The Weibull distribution can be generalized to obtain 
a Weibull regression model by allowing the failure rate to be 
a funeeton Of a vector Of qualitative personnel 
characteristics, X. In this application X is a row vector 
consisting of the binary variables X,j:0,, Xadaministratives Xtechnical# Xakilled # 
and Xecrfesion Which take On a value Of J Jt Sine observarided 
falls into the subscripted category, and are 0 otherwise. 


If the conditional (on X) hazard is 


X(t] X) =Ap (At)? te? 
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then the conditional density of T is 


f(t|X) =Ap(At)? tee tare’) (4.4) 


where 6=(6,,...,6,) is a vector of regression parameters. In 


terms of Y=log T, equation (4.4) is the linear model 


Y=a + XB*+ oW (4.5) 


where a=-log A, o=p', B'=-08, and W has the extreme value 
distribution. In equation (4.5) the personnel characteristics 
act additively on the time scale, vice multiplicatively in 
equation (4.4). [Ref. 1] 
3. Modelled Survivor Function 
Equation (4.5) is known as an accelerated failure time 


model. It is characterized by the scaling of the failure time 


Laer eT 


where T, 1s a failure time from a baseline when the values of 
X are all zero. Using this time scaling the modelled survivor 


function is 


Z5 


S(t) =e er’ (4.6). 


In equation (4.6), S(t) can be estimated by inserting 
appropriate values for X and t, and estimated values of 8B, \X, 
and p. This will produce a modelled survivor function which 
can be plotted to produce curves like those in Chapter III. 
If this modelled survivor function is similar (in shape and 
magnitude) to the empirical survivor function, then the model 
can be used. If it is not similar, then it possibly should 


not be used. 


C. RESULTS 

The partitioned data sets were regressed using the SAS 
LIFEREG procedure [Ref. 4] to obtain estimates for 8, r, aaa 
p. These values along with appropriate values for X were 
inserted into equation (4.6) to estimate modelled survivor 
funcerens: Figures 8 through 13 on pages 43 through 46 
contain plots of these modelled functions against empirical 
Survivor functions. Generally the modelled survivor functions 
do not fit the data as well as desired. 

Nonetheless, modelled survivor functions for these two 
paygrades will be used pending the development of sharper 
alternatives. The partitioned data sets are so small, it is 
felt that a model will produce the most reliable and stable 


probability estimates. 
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V. BINOMIAL COUNTING MODEL 


A. INTRODUCTION 

The Binomial Distribution provides a means for calculating 
expected values and variances of sums of random variables that 
have one of two outcomes: success or failure. This study is 
concerned with estimating the number of enlisted who separate 
on any given month (or conversely the number who do not 
separate). This chapter develops a simple model, based on 
principals of the Binomial Distribution, that calculates 
monthly survival figures per paygrade, and provides an example 
of how the model works. The primary reference is Chapter 3 of 


Reference 5. 


B. BINOMIAL DISTRIBUTION 
A binomial experiment is an experiment that possesses the 
following properties: 


- the experiment consists of n identical trials, 


- each trial results in one of two outcomes, a success or 


failure, 


- the probability of success on a single trial is equal to 


p and remains the same from trial to trial, 
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- trials are independent, and 


- the random variable of interest is Y, the number of 
successes in ntrials. The expected value of Y is 
E(Y)=np, where n is the number of trials and p is their 


corresponding probability of success. [Ref. 5] 


A person separating from the CG meets these standards and 
can be considered a Bernoulli trial. Each term trial results 
in one of two outcomes. A success is marked if the person 
does not separate, a failure when the person leaves the CG. 
The product-limit-estimate of the survivor function provides 
the probability of a person surviving to and beyond a time, t; 
P(Tet) for t=1,...0. This probability estimate remains the 
same for each trial (ie. for all persons with length of 
services equal to time t). The decision by a person to 
separate (or the CGs decision to separate a person) is not 
affected by another persons choice and thus independent of any 
other persons decision. And finally, Y equals the number of 
persons not separating within a paygrade who have survived 
beyond a time t. Using these definitions a basic model can be 
constructed that generates an expected value and variance for 


the number of persons that do not separate on a given month. 
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C. COUNTING MODEL 

Paygrades are partitioned in two steps to form data cells 
that consist of independent and identical trials required of 
the Binomial Distribution. The first step separates each 
paygrade into the single personnel characteristic examined in 
the survival analysis. These will be called sets. The second 
step further partitions the paygrade by the amount of time 
each person has served, or Time in Service (TIS) cells. 
Figure 14 provides a tree diagram example of this partitioning 


scheme. 


PAYGRADE L- 


CELL2 CELLr 


aee 
CELE1 CELL2 CELLr 


Figure 14. Counting Model Partitioning Scheme 





T is the maximum TIS any one person has within each set. 


Paygrades E-1 and E-2 are partitioned by race to form two data 


Zo 


sets each; paygrade E-3 is partitioned by MOS obtainment; and 
paygrades E-4 through E-9 are partitioned by MOS group. Each 
cell has a corresponding probability (ef @success isp) a curicr 
estimated from the actual data or generated by a model) which 
was defined in Section 4.D.1 to be the probability of survival 
being greater than or equal to t. Given@ay Cie wexpeered 
number of persons who fall into cell j, the expected number of 


successes (persons that do not separate) in each cell is 
E(Y) =p, 7 912 eee (a7) 
The cell variance is 


Var (Y,) =n,p,q, (4.8) 


where qj=1-p,;. The expected number of successes within each 


set is the sum of the cell expected values 


T 
E(Y,,) => E(Y,) (4.9) 
j=l 


a0 


and by summing cell variances the set variance is 


ie 
Var (Y,,) =}> Var(¥) (4.10) 


The expected number of successes per paygrade is calculated by 
summing set expected values. Paygrade variance is the sum of 
all set variances. 

Appendices C and D contain procedures for generating 


product limit estimates which correspond to the pjs. 


D. EXAMPLE 

Data from paygrade E-1 is used to illustrate this counting 
model. Paygrade E-1 contains 2 paygrade sets. Set, contains 
records for caucasians and Set, contains records for non- 
caucaSians. Table 6 and Table 7 on page 32 contain frequency 
counts of the number of persons falling within each cell per 


See 
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TABLE 6. PAYGRADE CELL FREQUENCY COUNT CAUCASIAN E-1 E- i 


fence | Pregame) Seok Freq 





TABLE 7. PAYGRADE CELL FREQUENCY COUNT NON-CAUCASIAN E-1 -CAUCASIAN E- “i 





Tables 8 and Table 9 on page 33 contain the eStimated p; values 


for each data sets. 


a2 


TABLE 8. PRODUCT LIMIT ESTIMATES PAYGRADE E-1 CAUCASIAN 

| vontn | op | month |p | month | oP 
aie ore per a Mawes | ss | 
a ee a ey ee ee ee 
Pere ree oe esse Yas a2! 
SS Ee ee a ee eee 
ee ee 0 ea ee eee 
ae ee ae a ee 
~417 -294 
a Ss a eo 
ee ee ee ee eee 















N |N 


THT 


ja 









TABLE 9. PRODUCT LIMIT ESTIMATES PAYGRADE E-1 NON-CAUCASIAN 





oS 


Substituting n; the values from tables 6 and 7 for the n, 
values and p, the estimated probability values from Tables 8 
and 9 into equation (4.9), the estimated expected value of 


each set is 
m 60 
E(Y,)=}> A,6,=353.475 
j=l 


52 
E(Y,) =)> Ap,=76.905 
j=l 


where set, corresponds to caucasians and set, corresponds to 
non-caucasians. Using equation (4.10) the Var(Y,) = 129.5 and 
Var(Y,) = 30.97. The overall paygrade estimated expected 


value and variance is the sum of the set values 


2 
E(Y,.,) =)> E(¥,) =430.38 
i=1 


2 
Var (¥,.,) =). Var (¥,) =160.47. 


i=1 
For working purposes these values are rounded down. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

There are approximately 30,100 enlisted persons in the CG 
at any given time. A method is needed to predict monthly 
attrition rates to maintain the CG enlisted hierarchical 
structure. Survival analysis was used to examine the 
characteristics of sex, marital status, race, and MOS group 
behavior of CG enlisted personnel. The results were fairly 
interesting. Males and married persons tended to have higher 
probabilities of survival than did females and individuals who 
were not married. Non-caucasians in most of the paygrades had 
higher survival probabilities than did caucasians. MOS groups 
demonstrated different characteristic behaviors among 
different paygrades. MOS groups in paygrades E-4 and E-5 
behaved about the same up to 48 months. After 48 months all 
groups showed a rapid decline in survivor probabilities with 
the aviation group having the highest probabilities and the 
technical group having the lowest. This drop is indicative of 
the end of a four year enlistment when individuals can usually 
separate from the CG. MOS groups in paygrades E-6 through E-9 
all had very high probabilities of survival and demonstrated 
a decline at certain months. E-6 showed this decline at about 


200 months and E-7 through E-9 demonstrated this at about 240 
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months. The drop in probabilities at 240 months is easily 
explained. Many Coast Guardsmen retire at 20 years and 
receive a substantial pension and many important benefits (ie. 
medical coverage for the person and his or her dependents). 
The drop at 200 months in paygrade E-7 is an anomoly. An E-7 
at this point in his/her career has almost 17 years of service 
time and has only about three more years to serve before 
retirement. Separating at 17 years seems to be financially 
inappropriate because the person will not receive a pension 
nor benefits if they separate. 

Modelled survivor functions were developed for paygrades 
E-8 and E-9 because these paygrades possessed small population 
sizes. Despite the fact that the models produced generally 
poor fits, it was felt that these survivor functions would be 
more reliable than survivor functions generated from very 
small data sizes. Therefore models were used to calculate 
attrition in paygrades E-8 and E-9. 

The binomial counting model provides an effective tool for 
calculating monthly attrition figures. Expected values and 
variances can be calculated by partitioning each paygrade into 
several binomial experiments. This model is easy to use and 


can be expanded in future studies. 


B. RECOMMENDATIONS 
As a first study the goals of this thesis were met. 


However further studies are required to better analyze CG 
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enlisted attrition. The counting model could be adjusted or 
refined based on findings from the recommendations listed in 
this section. 
1. Recommendation 1 
The counting model needs to be verified. Personnel 
data from early months of FY 92 need to be processed through 
the model to verify that it provides better results than the 
current calculation method. 
2. Recommendation 2 
Additional historical data needs to be analyzed to 
provide temporal stability to the study. The counting model 
produced in this study assumes that the characteristic 
behaviors of individuals identified in FY 91 will repeat 
themselves in FY 92. 
3. Recommendation 3 
A more manageable data base needs to be identified 
that allows for accessible man-machine interface. An 
accessible data base allows the analyst to verify, manipulate, 
and easily examine data entries aS an on-going task. This 
provides data reliability. The CG PPC is not responsible for 
Maintaining this type of data base, thus PWP should examine 
this issue. 
4. Recommendation 4 
Individual MOSs within MOS groups need to be examined 


to ensure they demonstrate similar attrition behaviors. 


a4) 


5. Recommendation 5 
Combined personnel characteristics need to be analyzed 
for survival behaviors and seasonal attrition trends need to 
be examined. 
6. Recommendation 6 
The data base needs to be extended in time. At 
present there is no justification for assuming time 
homogeneity of rates. Such an extension would also allow the 
study of some important exogenous variables, such as 
unemployment rate and inflation. 
7. Recommendation 7 
The educational levels of Coast Guardsmen should be 
analyzed for their effect on attrition behaviors. High school 
graduates, alternate high school degree holders, and non-high 
school graduates need to be considered to determine if there 


are differences in attrition behaviors amongst these groups. 
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APPENDIX A. PROCEDURES FOR ACCESSING DATA 
Data is stored on two magnetic disks. To access disks: 


a. login to mainframe computer 
b. type "getmvs mvs215 (separation data tape)" or 
"getmvs mvs217" (active duty data tape) 


To run a program that has a huge output file, ie. a file 
that will not fit on an A-disk: 


a. type "tdisk 6" to get a temporary disk with 6 
cylinders 

b. switch A-disk with temporary disk ( say temp disk was 
assigned "D 130" 
1. type "access 130 A" 
Zee ype access 191°)" 

Cc. copy program from the now D-disk to the newly created 
A-disk 

Ce LUM o@og ram 

e. switch disks back when finished by reversing procedure 


vee 


Be 


APPENDIX B. FIGURES 
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Figure 1. Survivor Function Paygrade E-1 by Marital Status 
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Figure 2. Survivor Function Paygrade E-4 by Sex 
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PAYGRADE E-1 BY RACE 
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Figure 3. Survivor Function Paygrade E-1 by Race 
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Figure 4. Survivor Function Paygrade E-2 by Race 
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SURVIVAL FUNCTION PAYGRADE E-3 BY RATE 
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Figure 5. Survivor Function Paygrade E-3 by Rate Obtainment 


SURVIVAL FUNCTIONS BY RATE GROUP 





Figure 6. Paygrades E-4, E-5, E-6 by Rate Group 
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SURVIVAL FUNCTIONS BY RATE GROUP 
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Figure 7. Paygrades E-7, E-8, E-9 by Rate Groups 
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Figure 8. Model Fit Paygrade E-8 
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Figure 9. Model Fit Paygrade E-8 
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Figure 10. Model Fit Paygrade E-8 
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Figure 11. Model Fit Paygrade E-9 
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Figure 12. Model Fit Paygrade E-9 
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Figure 13. Model Fit Paygrade E-9 


46 


APPENDIX C. PROCEDURES FOR GENERATING SURVIVAL FUNCTIONS 


FILE: LIFETEST CODE2 Al 


TITLE*E-2 BY RACE’ /e® THIS SAS CODE TAKES A DATA SET AND GENERATES 
SURVIVOR FUNCTIONS FOR THE DATA. THIS CODE WAS USED FOR PAYGRADE E-2 
AND CREATES SURVIVAL FUNCTIONS FOR CAUCASIANS AND NON-CAUCASIANS. 9/3 


DATA P02) 


INPUT RATE RANK SEX AGE MAR RACE MONTHS CEN /®THIS LINE DEFINES THE 
VARIABLES WITHIN THE DATA SET.®8/3 


GROUP s(RACE*S) /*THIS SEPARATES ALL CAUCASIANS FROM NON-CAUCASIANS. 
$ INDICATES CAUCASIAN. THIS SPLITS THE DATA SET INTO 2 SET? I CONTIANS 
CAUCASIANS. THE OTMER CONTAINS ALL OTHER MINORITITES. 8/3 


CARDS /*TELLS SAS TMAT THE DATA SET FOLLOWS. DATA SET HUST BE PLACED 
PHREOIATELY AFTER THIS LINE AFTER THE SEMICOLON, AND MUST HAVE A 
SEMICOLON PLACED IMMEDIATELY AFTER THE LAST FILE.®/ ; 


i 
PROC LIFETEST PLOT#(S) MAXTIME S48 OUTSERV*POCEST /*TMIS IS THE MAIN CODE 
IT CREATES A SURVIVAL FUNCTION PLOT (S). AND LIMITS THE X-AXIS OF THE 
PLOT TO <8 MONTHS. ALSO STGRE THE PRODUCT LIMIT ESTIMATES IN A S43 FILE 
CALLED PCCEST.«/ ; 


TIME MOMTMS*CEN(1) /®@DEFINES THE TIHE VARIABLE (MCNTHS). AND SETS THE 
CENSCRING VARIASLE NITHIN THE DATA SET EQUAL TO 1 ®/ 4 


STRATA GROUP /*@PERFORMS THE PROCEDURE ON ALL DEFINED GROUPS. IN TMIS CASE 
TRERE ARE ONLY 2 GROUPS®/ i 


PROC PRINT DATAsPOCEST /@PRINTS THE DATA SET POZEST INTO A CHS FILE®/: 


FILE: LIFETEST CODE al 


TITLE'PAYGRADE E~9 BY RATE GROUP ° /@SAME AS PREVIOUS CODE EXCEPT 
IT DIVIDES THE PAYGRADE INTO $ GROUPS: 1 FOR EACH MOS GROUP. GROUP I 
15 THE AVIATION PEOPLE, 2*ADMIN »SeTECHNICAL. “sSKILLED.5=CG PRO @8/ 3 


DATA PO9: 
INPUT RATE RANK SEX AGE MAR RACE MONTHS CEN: 
IF RATE=S60IRATE=S70IRATE2S20IRATESSSOIRATE*SS0 THEN GP#li 
ELSE IF RATE=3601RATE #4620 THEN GPt2i 
EU CE slr RATE S7S90IRATESIBOIRATESISOIRATE2270 IRATE S240 1RATE*280 THEN GPs33 
ELSE 1F RATE=SOOIRATE:84601RATE2870 THEN GPs; 
E\.cE GPsS: l 
CARDS: 
i 
PROC LIFETEST PLOT#($) MAXTIHE®sSO0 OUTSERV=POSEST: 
TIME HONTKS*CENCI D3 
STRATA GP: 
PROC PRINT DATASPOSESTS \ 
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eNO gna nannnnnnnnnnnnnaa 


INTEGER I.M.NLP 
PROGRAM MODFIT? 
REAL LAMBDA] .P].BETA]. BETAS, BETAS. BETAG.PETAS BETA] S$ 
RE AL LAMBDAD.PO.BETAG.BEYA?. BETAS, BETAY.BETAIO.BETAII.E 


THIS PROGRAM GENERATES MODELLED FRODUCT LIMIT ESTIMATES FOR ' 
REAL SHAT] .SHATS.SHATS.SHATC.SHATS .LAMBDAS,P3 


PAYGRADES WHICH MAVE BEEN PARTITIONED BY $ MOS GROUPS: AVIATION. 
ADMINISTRATIVE. TECHNICAL, SKILLED. CG PROFESSIONAL. (T USES AS 

INPUT REGRESSION COEFFICIENTS (BETA] THROUGH BETAIS). SHAPE 
PARAMETERS (P] AND P2), SCALE PARAMETERS (LAMBDAI AND LAMBDA2), 

AND A TIME PARAMETER (T). THIS PROGRAM IS ALSO DESIGNED TO WORK PARAMETER (MeIO-LAMBDAI®.1007 .P102.6B5) 

ON DATA THAT FIT A 3 PIECE WEIBULL. VARIABLES BETAI THRU BETAS, PARAMETER (BETA1©0.0.BETAZ0-.1938 .BETASe-. 0833, BETAG2-, 8 
Pl. AND LAMBDAI ARE VALUES FOR THE FIRST WEIBULL PIECE (P IS THE PARAMETER(N®200 -LAMPDAD?0.00003 .Pl06.638) 

TIME VALUE OF WHERE THE FIRST BREAK OCCURS). VARIABLES BETAS THRU ' PARAME TERI BETAS ©-.05<65 .BETAGe.119.BETA70~, 0519. BETAR e~< 
BETAS. P2. AND LAMBDA2 (AND M) CORRESPOND TO THE SECOND PIECE. PARAMETER (PeG20-LAMEDAS®.0028 .PS*8.0336 ) 

VARIABLES BETAS THRU BETAI2. P3. LAMBDAS (AND N) CORRESPOND TO THE PARAMETER(BETA9®.0256.BETAION. 0307. BETAI1°.106. BETA] 20- 


c FOLLOWING ESTABLISHES ALL VARIABLE VALUES. 


LAST WEIPULL PIECE. VARIABLE BETAI3 EQUALS ZERO AND CORRESPONDS PARAMETERIBETAIS20.0 ) 
TO THE CG PROFESSIONAL COEFFICIENT. THE VARIABLE SHAT I$ THE 

MODELLED PRODUCT LIMIT ESTIMATE: SHAT] FOR THE AVIATION MOS C CALCULATES SHATS FOR FIRST WEIBULL PIECE. 
GROUP, SHAT2 FOR ADMIN GROUP, SHATS FOR TECHNICAL GROUP. SHATS FOR 

THE SKILLED GROUP. AND SHATS FOR THE CG PRO GROUP. SHAT VALUES AR DO 10 1el.m 

CALCULATED FOR TIME PERIODST. (Tels...+ N)+ WHERE N IA THE LAST - SHATI ©EXP (= ((LAHBDAI © eEXP(-BETAL))@°0P1)) 


TIME VALUE NEEDED FOR THE PAYGRADE. a SHATS SEXP (-( (LAMBDA] ©] SEXP (-BETA2) )aePy])) 
SHATS*EXP (- ( (LAMBDA ® [EXP (-BETAS) )eaP})) 


THE ESTIMATES ARE WRITTEN ON THE CMS FILES: FISSFOO!. FYOSFOO!. 
SHATGCEXP (-((LAMBDAI ©] EXP (-BETAG ))aeP] )) 


AND FTSO0FOO]. 
SHATS EXP (~ ((LAMBDAI ©] EXP (-BETAIS) )eeP])) 
WRITE (886,91 )1. SHAT] »SHAT2,SHATS .SHATS.SHATS 
10 CONTINUE 


c CALCULATES SHATS FOR SECOND WEIBULL PIECE. 


DO 20 I2:M*].N 
SHATISEXP (~ ((LAMPDAL SI SEXP (-BETAS ) )aaP2)) 
SHATO SEXP (-( TLAMBDALO] SEXP (-BETAG) )eeP2)) 
SHATS SEXP (= ((LAMBDA2 ©] *EXP(~PETA?) )eePs)) 
SHATGZEXP(~( (LAMBDAL ©] SEXP (-BETAS) JeeP2)) 
SHATS 2EX> (- ( (LAMBCAL®] SEXP (-BETAIS) )eeP2)) 
WRITE (89.911, SHAT] -SHAT2.SHATS. SHATC. SHATS 
20 CONTINUE 


c CALCULATES SHATS FOR LAST PIECE. 


DO 30 leNe],P 
SHATI SEXP(- ((LAMBDAS®] SEXP (~BE TAS) )aaP3)) 
SHAT2 SEXP (-( (LAMBDAS SI SEXP (-BETAIO) )aePs)) 
SHATS SEXP (~( (LAMBDAS®] SEXP (-BETAII) )aePs)) 
SHATG CEXP (~ ( (LAMBDAS®] SEXP (-BETAI2) )eePs)) 
SHATS EXP (-( (LAMBDAS ©] CEXP(~BETAIS) )eePS)) 
WRITE(9O0,91)1.SHATI. SHAT2.SHATS.SHATC, SHATS 
30 CONTINUE 
$I FORMATCIX. IS. SK .FO GSK eFO Se SKeF GE. 6.3K LFE 63K FEC) 
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APPENDIX D. PROCEDURES FOR CONDUCTING REGRESSION ANALYSIS 


PROGRAM REGPREP 


Cc THIS CODE SEPARATES A PAYGRADE DATA SET BY MOS GROUP AND CREATES Cc THIS 1S THE LAST PIECE OF THE WEIBULL. 
Cc DATA SETS NEEOED TO CONDUCT THE SAS LIFEREG PROCEDURE. AGAIN THIS 
c CODE 1S DESIGNED FOR A 3 PIECE WEIBULL. THE INPUT FILE (#2) ELSE 
Cc CONTAINS THE PAYGRADE DATA SET. 3 OUTPUT DATA SETS ARE CREATED TO IF(RATE .EQO. $20 .OR. RATE .EO. $60 .OR. RATE .£0.870 .OR. 
Cc CONFORM TO THE § PIECE WEIBULL. THE VARIABLE N 1S THE NUMBER OF 1 RATE .EQ. $30 .OR. RATE .EQ. SSO) THEN 
c FILES 1N THE PAYGRADE. WRITE(S8.93) HONTHS.CEN 
ELSEMF(RATE .E9. 360 .OR. RATE .EQ. 620) THEN 
INTEGER 1.N.RATE. MONTHS CEN WRITE([8.96) MONTMS.CEN 
PLRAMETER(N2668) ELSEIF(RATE .£0. 270 .CR. RATE .EQ. 260 .OR. RATE .EO. 180 
OFEN(2.FILE®'POS DATA A’) ; 1 OR. RATE .EQ. 790 .OR. RATE .EQO. 1460 .OR. RATE .EO. 
DO 10 1*).N . 1 280) THEN 
READ(2,92) RATE.HONTHS.CEN : WRITE (28.95) MONTMS,CEN 
ELSEIFC(RATE .EQ. 870 .OR. RATE .E0. $00 .OR. RATE .E0. 840)THEN 
c THIS 18 THE FIRST PIECE OF THE WEIBULL WRITE (28.96) MONTHS.CEN 
ELSE 
1F (MONTHS .LE. 10) THEN WRITE (28.97) MONTMS.CEN 
IF(RATE .EQ. $20 .OR. RATE .£0. S60 .OR. RATE .£0.870 .OR, ENDIF 
} RATE .EQ. 530 .OR. RATE .EQ. SSO)THEN ENDIF 
WRITE (26.93) MCNTHS.CEN 10 CONTINUE 
ELSEIF(RATE .EQ. 360 .OR. RATE .£Q0, 420) THEN 
WRITE (26.94) MONTHS.CEN Cc THIS 1S THE FORMAT FOR THE READ/INPUT STATEMENTS. 
ELSEIF(RATE .£0. 270 .OR. RATE .EQ. 260 .OR. RATE .E0. 180 
} -OR, RATE .EQ. 790 .OR. RATE .EQ. 140 .OR. RATE .EO. 
: ee } 92 FORMATCIX,DSe13X+1561K017) 
WRITE(26.9S) MONTHS,.CEN 
ELSEIF(RATE .E0. @70 .OR. RATE .E2. S00 .OR. RATE .£0. ®60)THEN - THESE FORHAT STATEMENTS CREATE REGRESSION DATA SETS. MONTHS 1s THE 
WRITE(26.96) MONTHS CEN . DEPENDENT VARIABLE AND 1S LOCATED IN THE FIRST COLUMN, THE 
Poe ¢ PRECEOING COLUMNS CONTAIN A 0 OR I IF THE PERSON FALLS INTO oe 
aa ) €——sCOTF THE MOS GROUPS. (THE STH MOS GROUP 1S ASSUMEO AND DCES NO 
ee C IN THE DATA SET. 
c THIS 1$ THE SECOND PIECE OF THE WEIBULL. 93 FORMATCIX. 130 1Ke Do 1K OT THs OTe DKe Oo IMeDT) 


} 96 FORMATCINs 1301 OTe IX Le IKe Oo IMO IKI) 


ce Oe Qe ee > EE, 200) THEN 9S FORMAT(IMs 1301007 DMs Ooo IM e Eo DMe OTK TT) 


IFIRATE .E0. S20 .OR. RATE .EQ. S60 .OR. RATE .EO.570 .OR. 


’ 96 FORMATCIX, TS «160% 1K 6 OT TKS OTe IMA TAMIL) 
! RATE .£0. S30 .OR. RATE .£0. SSO) THEN ee cia water ix POTN. 10 o1Nes @ys1N- 10) 
WRITE(27.93) MONTHS CEN at 
ELSEIF(RATE .£0. 360 .OR. RATE .£0. 620) THEN Sra 
WRITE(27.94) MONTHS «CEN 
ELSEIF(RATE .£0. 270 .OR. RATE .£0. 240 .OR. RATE .£0. 180 
,OR. RATE .EQ. 790 .OR. RATE .£0. [<p .OR. RATE .€0, 
I 280) THEN 


WRITE(27.95) MONTHS.CEN 
ELSEIF(RATE .EQ. 870 .OR. RATE .EQ. S00 .OR. RATE .EQ. 840) THEN 
WRITE(27.96) MONTHS.CEN 
ELSE 
WRITE(27.97) MONTHS CEN 
ENOIF 
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' 


FILE: LIFEREG CODE Al 


TITLE REGRESSION OUTPUT PAYGRADE E-6' /=SaS CODE PERFORMS A PEGRESSION 
ON CENSORED DATA. GENERATES REGRESSION COEFFICIENTS. SHAPE AND SCALE 
PARAMETERS, OUTPUT 15 VERY SMALL. ®/ 3 


DATA POG; 

INPUT MONTHS AVIATION ADMIN TECH SKILLED CEN: 
CARDS: 

' 

PROC LIFEREG OUTESTsMODEL: 


CLASS AVIATION ADMIN TECH SKILLED /®SETS THE VARIABLES AVIATION, ADMIN, 
TECH. AND SKILLED AS CATEGORICAL (0 OR 1) VARIABLES */ 13 


MODEL MONTHS®CEN(] © AVIATION ADMIN TECH SKILLED/DIST#WEIBULL /*SETS THE 
DEFENDENT VARIABLE EQUAL TO MONTH, SETS THE CENSORING VARIAPBLE-CEN®]. 
SETS THE INDEPENDENT VARIABLES EQUAL TO AVIATION. ADMIN, TECH. SKILLED. 
ESTABLISHES THE WEIBULL AS THE PARAMETRIC DISTRIBUTION ®/ 5 
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